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We investigate magnetic order in a lattice of classical spins coupled to an isotropic gas of one¬ 
dimensional (Id) conduction electrons via local exchange interactions. The frequently discussed 
Ruderman-Kittel-Kasuya-Yosida (RKKY) effective exchange model for this system predicts that 
spiral order is always preferred. Here we consider the problem nonperturbatively, and find that such 
order vanishes above a critical value of the exchange coupling that depends strongly on the lattice 
spacing. The critical coupling tends to zero as the lattice spacing becomes commensurate with the 
Fermi wave vector, signalling the breakdown of the perturbative RKKY picture, and spiral order, 
even at weak coupling. We provide the exact phase diagram for arbitrary exchange coupling and 
lattice spacing, and discuss its stability. Our results shed new light on the problem of utilizing a 
spiral spin-lattice state to drive a one-dimensional superconductor into a topological phase. 

PACS numbers: 75.30.Hx, 75.75.-c, 71.10.Pm, 03.67.Lx 


Introduction — The RKKY coupling mechanism serves 
as a cornerstone of our understanding of the indirect 
exchange interaction between magnetic impurities in a 
metallic host [1-3] . Recently, the RKKY effect has played 
a central role in the intense effort to realize topological 
superconductivity and associated Majorana bound states 
in systems of magnetic impurities placed on s-wave super¬ 
conductors [6-12]. In a similar setup, promising evidence 
for Majorana states was recently reported [17, 18] from 
experiments in which a scanning tunneling microscope 
was used to create and probe a chain of ferromagnetically 
aligned Fe atoms on the surface of bulk superconducting 
Pb. Although the physics of this specific realization is 
likely dominated by direct Fe-Fe exchange interactions, 
spin-orbit coupling and magnetic anisotropy, these excit¬ 
ing results provide motivation to develop a deeper under¬ 
standing of possible “self-organized” phases of magnetic 
impurities embedded in electronic environments. 

For a one-dimensional (Id) system without spin-orbit 
coupling, the emergence of topological superconductiv¬ 
ity is directly tied to the presence of spiral order in the 
magnetic impurity chain. Such ordering has been pre¬ 
dicted to arise naturally in the system’s ground state 
due to the RKKY coupling between magnetic impuri¬ 
ties [6-9]. Focusing on the normal state of the host, 
the RKKY coupling arises at second order in the di¬ 
rect exchange between conduction electrons and impu¬ 
rities. Because the electron wavefunctions oscillate in 
space with the Fermi wavevector kp, the mediated inter¬ 
action exhibits characteristic 2kp oscillations: FIrkky cx 
— Si • S 2 cos( 2 fcpri 2 )/ri 2 , where Si 2 represent the spin 
vectors of magnetic impurities separated by a distance 
ri 2 . The long-range ~ 1/r scaling of the interaction 
results from the gapless nature of excitations near the 
Fermi points of the Id host system, and is responsible 
for the logarithmically diverging static spin susceptibil¬ 
ity x(fc) (Fourier transform of FIrkky) at wavevector 
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FIG. 1. (Color online) Electronic band structure for the spin 
lattice in the antiferromagnetic state (left panel) and spiral 
phase (right panel). The corresponding points in the phase 
diagram are shown by the symbols 4|k and ^ in Fig. 2. Close 
to the commensurate points, kpa = (2n-|-l)7r/2, the AF state 
fully gaps the electronic system. When the lattice spacing is 
increased the gaps in the AF state become smaller and at a 
critical value of the lattice spacing the partially gapped spiral 
state becomes the lowest energy state. 

k = 2kp. 

Considering a classical spin lattice subject to this in¬ 
teraction, one easily sees that it lowers its energy E by 
forming a spiral state with • Sj = cos{qrij): 

E{q) oc - E cos{qrij)cos(2kprij) / r^j 

iAl 

=-Rex(g) ~ log|(2fcF - kl)a|, (1) 

where a is the lattice spacing. Minimizing Eq. (1) gives 
q = 2kF, with a logarithmically diverging energy gain [4]. 

This diverging energy gain is an artifact of the Born 
approximation, where renormalization of electron propa¬ 
gation due to spin-lattice scattering is neglected. In par¬ 
ticular, such scattering exhibits resonant enhancement 
due to Bragg reflection when 2fci?a ~ irn for integer n 
and is associated with the opening of electronic band 
gaps at the Fermi energy. At the odd commensurate 
points 2kFa = (2n -I- l)7r, the system is an antiferro- 
magnet (AF) with a fully gapped electronic bandstruc- 
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FIG. 2. (Color online) Groundstate phase diagram of the 
model in Eq. (2) as a function of lattice spacing a and ex¬ 
change coupling J. Shading schematically represents the wave 
vector q of the spiral, lying between g = 0 in the ferromag¬ 
netic (F) phase and q = -k j a va the anti-ferromagnetic (AF) 
phase. Solid phase boundary lines correspond to first-order 
transitions, double lines to second-order transitions and dots 
to triple points. The symbols 4|k and correspond to the 
parameters used for the bandstructures shown in Fig. 1. 

ture, while a spiral phase (S) gives a single band gap, see 
Fig. 1. For any finite coupling, these two bandstructures 
are not smoothly connected. Hence the AF phase may 
be expected to persist throughout extended intervals of 
lattice spacings, around the commensurate values. The 
competition between partially and fully gapped S and 
AF phases demonstrates the relevance of higher order ef¬ 
fective n-body spin interactions, neglected in the 2-body 
spin Hamiltonian I?rkky- 

In this Letter, we investigate the ground state of a clas¬ 
sical spin lattice coupled locally via spin exchange to an 
isotropic gas of Id conduction electrons. We find that 
spin lattice spiral order vanishes above a critical value 
of the exchange coupling, which depends strongly on the 
lattice spacing. In particular, for 2kFa ~ mr where n 
is an even (odd) integer, the second (first) order phase 
transition to the ferromagnetic (anti-ferromagnetic) state 
occurs even at small exchange coupling. Remarkably, our 
exact calculation shows that the cumulative effect of all 
n-body interactions (neglected in the RKKY/Born ap¬ 
proximation) can drive the system into new phases, even 
for weak coupling where the RKKY picture is naively 
expected to be valid, see Fig. 2. 

Similar competing magnetic phases were reported in 
Refs. [19-21], which focused on the dense impurity limit 
kpa < 1 in higher dimensions [19, 21] or in the long 
wavelength limit go <C 1 [20] . Motivated by the prospect 
of achieving a topologically non-trivial phase using mag¬ 
netic adatoms controllably arranged on a superconduct¬ 
ing surface, we consider a continuously variable spin- 
lattice spacing. This, in particular, allows us to address 
the dilute impurity limit, kpa > 1, and to show that 
commensurability plays an important role for arbitrary 
kpct. 

In addition to exposing the competing phases, our 


analysis allows us to determine the exact phase bound¬ 
aries at arbitrary exchange coupling J and lattice spac¬ 
ing. Interestingly, above a critical exchange coupling 
J > Jc ~ the spiral phase vanishes for any lat¬ 
tice spacing (here n = 2j {tt ^ is the chemical 

potential, and m is the effective mass of electrons in the 
metallic host). The existence of the critical coupling im¬ 
plies the appearance of triple points in the phase dia¬ 
gram, where all three phases coexist in thermodynamic 
equilibrium, see Fig. 2. 

Model — We illustrate the effects described above by 
studying the ground states of an electronic system with 
the Hamiltonian {h= 1) 


H = 


dx 


^ 2m 


+ Je^S{x) ■ S{x) , 


( 2 ) 


where ip a (V'J) space fermionic annihilation (cre¬ 

ation) operators with spin projection a, p = —idx, 
s = ^ipp^cr^^iip^i is the electron spin density (cr is the 
vector of Pauli matrices) and Jex is the exchange inter¬ 
action constant. We model the magnetic lattice by clas¬ 
sical spins, described by S(x) = ^{x — ja) Uj, where 

a is the spacing between impurity spins and rij is a unit 
vector denoting the direction of the spin. The quan¬ 
tum limit of a Hondo lattice of spins with ^ = 1 /2 has a 
separate, rich history [22-27]. 

To investigate magnetic ordering in the system de¬ 
scribed by Hamiltonian (2), we use a variational ap¬ 
proach in which we minimize the electronic ground state 
energy with respect to the magnetic ordering profile. 
We consider a planar spiral order ansatz of the form 
Uj = (cos qaj, sin qaj, 0) which interpolates between fer¬ 
romagnetic (F) {q = 0) and anti-ferromagnetic (AF) 
{q = ±n/a) order [31]. 

To evaluate the groundstate energy of Hamiltonian (2), 
we first remove the winding of the ordering field using 
the unitary transformation U = [5-9, 28]. The 

transformed Hamiltonian H' = UHW has discrete trans¬ 
lational symmetry, as well as an effective spin-orbit cou¬ 
pling in the direction perpendicular to the spiral plane, 



E 


{p- -ip a 

2m 


+ JexS'(a;) • s(a;) 


(3) 

where S' = SJ2j ~ 3^)^ is related to S by an SO(3) 
rotation around the 0 -axis. The Hamiltonian H' is rem- 
iniscient of models of quantum wires with intrinsic spin- 
orbit coupling in a homogeneous perpendicular magnetic 
field [13-16]. However, the discrete and self-ordering na¬ 
ture of the effective magnetic field produced by the spin- 
lattice are essential in our work, and lead to interesting 
new phenomena. 

Within the variational approach, we determine the op¬ 
timal wave vector q of the spin lattice by minimizing the 
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FIG. 3. (Color online) Optimal wave vector of the spin spiral 
q as a function of exchange coupling at kpa = O.Qtt (red) and 
kpa = O.dvr (blue), shown as vertical cuts of the inset. 


FIG. 4. (Color online) Optimal wavevector q and energy per 
impurity spin Eq = Eoa/L as a function of kpa at Jv = 0.2, 
shown as horizontal cut of the inset in Fig. 3. The energy is 
measured relative to the q — 0 state. 


groundstate free energy Eo{q) = {H — fiN)o for fixed /r 
{N is the total electron number). Up to a q-independent 
constant, E^i^q) can be written in terms of the shift in 
the electronic density of states, dv, due to the presence of 
the spin-lattice potential: i?o(9) = “ 

eSiy{e + fi). Using Lloyd’s formula [32, 33], we 
express the shift in the density of states in terms of the 
free electron Green’s function Gq and the spin-dependent 
potential V: Sv(e) = — ilmTr9£ln(l — Go(e)U). 

In order to evaluate the shift of the density of states, 
Siy(e), we exploit the discrete translational symmetry of 
H', which allows it to be block diagonalized using a set of 
Bloch states labeled by quasi momentum jfcj < n/a (mod¬ 
ulo ja). Importantly, the delta function lattice acts as 
a rank-1 operator within the subspace of states for each 
value of the quasi momentum k. Using the Woodbury 
matrix identity, we can thus reduce the operator trace 
occuring in 6v to an integral over the Brillouin zone (the 
remaining trace over the 2x2 spin subspace can be com¬ 
puted explicitly). In this way we obtain the ground state 
energy (per system length L) for each fixed q, measured 
relative to the value with J = 0: 


^o(g) 

L 



TT 



dk 


— oo _7r/a 



where Go{e,k) = = 

(TOa/a)sinaa [cos (fca =F ga/2) — cosa] ^ is the Green 
function for free electrons (i.e., evaluated for J = 0) 
with conserved spin projection ±1 along the z-axis, and 
a = sj2m{e -I- *0+ -I- /r). Minimization of Uo(q) gives 
q = goptimai which we denote as q for brevity. We asso¬ 
ciate q — 0 with F, q = ir/a with AF and intermediate 
values with S, see Fig. 2. 

As noted in Refs. [6-9], the physics underlying the or¬ 
dering of the spin-lattice can also be understood by ana¬ 
lyzing the corresponding electronic band structure. From 


that point of view, the advantage of the q = 2kp spiral 
comes from the opening of a gap for half the degrees of 
freedom at the Fermi surface, see Fig. 1, which lowers 
the energies of occupied states. This “spin-Peierls effect” 
is responsible for the existence of the stable spiral phase, 
but ultimately the influence of resonant Bragg reflection 
leads to a transition near points of commensurability. 

Phase diagram for Jv 1 — The pole structure of 
Eq. (4), and thus the electronic band structure, is deter¬ 
mined by the zeros of the argument of ln(l — ^GqQq). 
Substituting in the expressions for Qq gives 


1 = 


2 t2 o 

™ / sin^ aa 


(cos kpa — cos aa) (cos k^a — cos aa) ’ 


(5) 


where k± = fc ± q/2. In order for the spin-Peierls effect 
to yield the maximal energy gain, the gap which it opens 
should be centered at the Fermi energy: e = ±A. This 
relatively simple (heuristic) condition on the gap, when 
imposed on the exact spectral equation (5), provides an 
approximation for the wave-vector of the spiral order. 


ggapO = 2Re arccos -I- {'kJv/2)‘^ jcosfc^a]^ , (6) 


which agrees quite well with the exact numerical result 
of Eq. (4), see red solid lines in Fig. 4. 

Remarkably, the approximate relation in Eq. (6) cor¬ 
rectly shows that the RKKY result q = 2kp is recov¬ 
ered only in the limit J —>■ 0. More generally, Eq. (6) 
predicts a critical value Jc^'^ = jtanfciT’oj at which 
the optimal spiral wavevector goes to zero continuously, 
q ~ signalling a second order phase tran¬ 

sition to the F state, see Fig. 3. In the F state, the 
electronic states at the Fermi surface are spin-polarized. 
Although not exact, the phase boundaries predicted by 
Eq. (6) provide a good qualitative description of the F-S 
boundaries in Fig. 2. 
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The S state also competes with the AF order near half 
integer values of /cfo/tt. At these points, the system may 
become fully gapped, potentially lowering its energy even 
further. The transition roughly takes place when the AF 
band gaps close, although we generally find that when 
the AF state is the groundstate the electronic bands are 
always gapped: the system is a band insulator, see Fig. 1. 
Using Eqs. (5—6), this gives the Ist-order transition to 
the AF state at J > ^ |cotA:Fa| and = 

|arccos(7rJ f/ 2) « f (1 — Ji^)- Thus the jump of the 
order parameter is 6q/q = Jv, to leading order in Jv. 

We thus arrive at the following sketch of the phase di¬ 
agram for Jv <C 1. For J < ^min (|tanfcFfl|, \cotkpa\), 
the ground state wavevector is approximately given by 
Eq. (6); otherwise, q = 0 (F) or tt/o (AF). As a result, 
Jc becomes arbitrarily small when 2kpa ^ irn for integer 
n, implying that the validity of the RKKY approxima¬ 
tion is restricted to Jv/ \kpa — mr/2\ 1 (not Jv 1). 

This analysis also suggests the presence of triple points 
for J:/ ~ 1 which are verified by the exact results (Fig. 2). 

Phase diagram for Jv ^ 1 — In this case the phase 
diagram is dominated by the competition between F and 
AF configurations. In the limit kpa 1 the dominant 
contribution to the energy comes from occupation of the 
deep bound state bands formed by the strong spin-lattice 
potential. For kpa S> 1 the bound state bands play little 
role and the phase diagram is controlled by positive en¬ 
ergy bands of the scattering continuum. After a more in¬ 
volved analysis of the spectral equation (5), one finds that 
the F state exists only in narrow regions around points of 
(even) commensurability such that Jv < 1/ \kpa — 7m/2| 
for even integer n. 

Discussion — The structure of the phase diagram in 
Fig. 2 is stable against a variety of perturbations. Gen¬ 
erally, we find that modifications remain small as long 
as the associated energy scale is small compared to J/a. 
Perturbation theory can then be used to estimate changes 
in the phase boundary lines. 

For nonzero temperature T, thermal spin fluctuations 
tend to deplete the magnetic order, which can be stabi¬ 
lized at T > 0 either by a finite chain size or a magnetic 
field/anisotropy, see [6-9] for a detailed discussion on this 
point. Thermal fluctuations also modify the occupation 
of electronic states, which changes the phase boundaries 
in Fig. 2. The associated increase in free energy is mini¬ 
mized in the (band insulator) AF state, where the depen¬ 
dence is exponential SF(T) ~ The free energy 

increases faster for the gapless S and F states, going as 
v{Ep)T‘^. The exact density of states can be studied 
using Eq. (5) and shows that for kpa < 1 the F state 
is favored. Thus, the S region for kpa < 7r/2, Jv < I 
shrinks with growing T. 

The phase boundaries are also stable under inclusion of 
weak electron-electron interactions. In the presence of a 
local density-density interaction Hint oc g 


one can use first order perturbation theory to determine 
the shift of the phase boundaries. At kpa < 1 the F 
phase corresponds to a spin-polarized Fermi sea and is 
thus unaffected to leading order, while the energy in the 
S phase shifts by an amount proportional to g. The S- 
phase is thus diminished (enhanced) for g > 0 (g < 0). 
Determining the direction of the shift at the S-AF phase 
boundary requires further analysis. 

More interesting is the role of intrinsic SO coupling 
of the electron gas. By choosing the plane of the spiral 
to align with the SO axis [34], we see that the inclusion 
of intrinsic SO coupling merely shifts the ordering wave 
vector to a net value: q ^ q — qso, leaving the phase 
diagram unaltered [5]. Thus, if the parameters (J, a) of 
the system lie in the F phase of Fig. 2, the groundstate of 
the spin lattice will have the wave vector q — qsoy which 
acts to “screen” the intrinsic SO coupling and renders 
the electronic bands ferromagnetic in the spiral-rotated 
reference frame. Similarly, bands in the AF phase also 
have a well-defined spin projection in the same frame. 

This “spin-lattice screening effect” has important im¬ 
plications for the prospect of achieving a one-dimensional 
topological superconductor and associated Majorana 
bound states by a proximity induced s—wave pairing[6- 
12, 29, 30]. In particular, note that the topological phase 
can only occur in the S phase in Fig. 2 (if the pairing is 
local in space). This is because the spin rotation induced 
by the net ordering wave vector needs to be non-collinear 
at the two Fermi points, which is only the case in the S 
regions of Fig. 2. This complication may be overcome to 
some extent by the pinning effect of an easy-axis mag¬ 
netic anisotropy energy ^ — (iS^)^, which favors collinear 
order in the direction perpendicular to the substrate sur¬ 
face (an easy-plane merely pins the spiral plane). As 
such, easy-axis anisotropy competes with the tendency 
of the spin-lattice to form the planar screening spiral. 

The one-dimensional regime can be realized if fj,, T are 
smaller than the transverse level spacing of a quasi one¬ 
dimensional wire. Occupation of higher subbands leads 
to a dimensional crossover. Here we expect the S phase 
to be to suppressed, since the 2kp instability in magnetic 
susceptibility does not exist in higher dimensions for the 
normal state [4, 5]. In the superconducting state, how¬ 
ever, this suppression may not hold due to the presence of 
resonantly enhanced spin-exchange via Shiba states [35]. 
Understanding the role of such effects requires further 
investigation. 

Conclusion — The magnetic phase diagram of a spin 
lattice embedded in a one-dimensional conductor exhibits 
a rich structure that strikingly reveals physics beyond 
the RKKY/Born approximation. The simplified model 
(2) admits an exact solution, which allowed us to de¬ 
duce the precise phase boundaries of the spiral order and 
determine the limitations of the RKKY approach. Our 
results clarify the regimes where topological supercon¬ 
ductivity based on a “self-organized” spin-lattice state 


5 


can be realized. 

Acknowledgements — The research was supported by 
the Danish National Research Council, the Villum Kann 
Rasmussen Foundation, and by the People Programme 
(Marie Curie Actions) of the European Union’s Seventh 
Framework Programme (FP7/2007-2013) under REA 
grant agreement PIIF-GA-2013-627838. 


[1] M. A. Ruderman and C. Kittel, Phys. Rev. 96, 99 (1954). 

[2] T. Kasuya, Prog. Theor. Phys. 16 , 45 (1956). 

[3] K. Yosida, Phys. Rev. 106, 893 (1957). 

[4] Already for d = 2 the RKKY model instead predicts 
ferromagnetic {q = 0) or anti-ferromagnetic {q = rr/a) 
order [5]. This is a result of the slower ~ 1/r^ decay of 
the induced spin-spin exchange interaction and thus an 
absence of a divergent spin susceptibility a,t q — 2kp. 

[5] Y. Kim, M. Cheng, B. Bauer, R. M. Lutchyn and S. Das 
Sarma, Phys. Rev. B 90 , 060401(R) (2014). 

[6] B. Braunecker and P. Simon, Phys. Rev. Lett. Ill, 
147202 (2013). 

[7] J. Klinovaja, P. Stano, A. Yazdani and D. Loss, Phys. 
Rev. Lett. Ill, 186805 (2013). 

[8] M. M. Vazifeh and M. Franz Phys. Rev. Lett. Ill, 
206802 (2013). 

[9] 1. Reis, D. J. J. Marchand and M. Franz, Phys. Rev. B 
90 , 085124 (2014). 

[10] F. Pientka, L. 1. Glazman, F. von Oppen, Phys. Rev. B 

88, 155420 (2013). 

[11] F. Pientka, L. 1. Glazman, F. von Oppen, Phys. Rev. B 

89 , 180505(R) (2014). 

[12] Y. Peng, F. Pientka, L. 1. Glazman, F. von Oppen, 
arXiv:1412.0151. 

[13] R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. 
Lett. 105 , 077001 (2010). 

[14] Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 
105 , 177002 (2010). 

[15] J. Alicea, Rep. Prog. Phys. Usp. 75, 076501 (2012). 

[16] M. Leijnse and K. Flensberg, Semicond. Sci. Technol. 27 , 


124003 (2012) 

[17] S. Nadj-Perge, et. ah, Science 346, 602 (2014). 

[18] J. Li, et. ah, Phys. Rev. B 90, 235433 (2014). 

[19] M. Hamada and H. Shimahara, Phys. Rev. B 51 , 3027 
(1995). 

[20] D. Pekker, S. Mukhopadhyay, N. Trivedi, and P. M. Gold- 
bart, Phys. Rev. B 72 , 075118 (2005). 

[21] K. Pradhan and P. Majumdar, Europhys. Lett. 85, 37007 
(2009). 

[22] H. Tsunetsugu, M. Sigrist, and K. Ueda, Rev. Mod. Phys. 
69, 809 (1997). 

[23] G. Honner and M. Gulacsi, Phys. Rev. Lett. 78 , 2180 
(1997). 

[24] G. Honner and M. Gulacsi, Phys. Rev. B 58 , 2662 (1998). 

]25] N. Shibata and K. Ueda, J. Phys.: Condens. Matter 11 , 

R1 (1999). 

[26] D. J. Garcia, K. Hallberg, B. Alascio, and M. Avignon, 
Phys. Rev. Lett. 93, 177204 (2004). 

[27] S. A. Basylko, P. H. Lundow, and A. Rosengren, Phys. 
Rev. B 77 , 073103 (2008). 

[28] B. Braunecker, G. 1. Japaridze, J. Klinovaja, and D. Loss, 
Phys. Rev. B 82 , 045127 (2010). 

[29] M. Kjaergaard, K. Wolms, and K. Flensberg, 
Phys. Rev. B 85, 020503(R) (2012). 

[30] J. Klinovaja, P. Stano, and D. Loss, Physical Review 
Letters 109, 236801 (2012). 

[31] Other spin configurations with, e.g., an out of plane 
component or higher harmonics, always appear to be of 
higher energy, consistent with previous results [5, 9]. Note 
that due to the SU(2) symmetry of the problem, the ori¬ 
entation of the spiral plane is arbitrary. 

[32] P. Lloyd, Proc. Phys. Soc. Lond. 90, 207 (1967). 

[33] P. Lloyd and P. V. Smith, Adv. Phys. 21 , 69 (1972). 

[34] The orientation of the spiral is not fixed by the intrin¬ 
sic SO coupling since the latter can be eliminated via a 
spin-dependent gauge transformation, analogous to the 
relation between Eqs. (2), (3). In the rotated frame, the 
system exhibits rotational invariance and the plane of the 
spiral may be chosen arbitrarily. 

[35] N. Y. Yao, L. 1. Glazman, E. A. Dernier, M. D. Lukin, 
and J. D. Sau, Phys. Rev. Lett. 113, 087202 (2014). 



